Annual Report 


♦ 








/V l /3 ' &/ 


Under 

Contract # NAG -5 - 1010 
National Aeronautics and Space Administration 




Improvement and extension of a radar forest backscattering model 


Prepared by: 
David S. Simonett 
Yong Wang 


(AASA-CB-183259) I BBECVEBJ- A1 AMD EXI’EBSICB N89- 11292 

Ci A EA£AB F0BES1 £ ACKSCA11IJ 1KG ECBEL 
A tnud 1 Beport (California Bri?.) 38 p 

CSC! 02-F Unclas 

G3/43 0166487 


Department of Geography 
University of California Santa Barbara 


September 20, 1988 



ANNUAL REPORT 


1. Introduction: 

Research to-date has focused on modeling development and programming based on model 
components we proposed during the past several months and research progress made by Simonett’s 
team. In this report, we summarize the model components and programs (in C language under 
UNIX) finished to-date. These model components may help explain the contributions of various 
vegetation structural components to the attenuation and backscattering of vegetated surfaces to 
extract useful data concerning forest stands and their underlying surfaces for both the seawater-on 
and seawater-off cases. 

2. Model components*: 

There are six model components as shown in Figure 1, which are described as follows: 

1 Direct backscattering from ground surface (water-off case only) with attenuation from upper 
and lower canopies (short form dbs). 

2 Direct volume scattering from upper tree canopy (short form vsu ). 

3 Direct volume scattering from lower tree canopy with attenuation from the upper tree canopy 

(two way) (vsl). - •* 

4 Interaction of ground / trunk forward reflection with attenuation from upper and lower tree 
canopies (two way) (tg). 

5 Interaction of ground / lower tree canopy forward scattering with attenuation from upper tree 
canopy (two way) ( glc ). 

6 Interaction of Ground / upper tree canopy forward scattering with attenuation from lower tree 
canopy (two -way) {guc). 

3. Programs: 

head.h: a header file in which all the constants and external variables are defined. These constants 
and variables will be used in the whole program of modeling radar backscattering from forest 
(Sundri and Gewa) stands for the Bangladesh research project. 


♦note: Upper canopy is the canopy of trees whose dbhs are >= 5.0 cm. Lower canopy is the canopy of 
trees whose dbhs are < 5.0 cm, which are also called as regeneration trees. 
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main.c: the main program, detailed as follows: 

NAME 

SIGMA - radar returns of forest stands. 

SYNOPSIS 

SIGMA intab indata inpxl indbh sts_o t_o dbs_o vsu_o vsl_o tg_o glc_o guc_o 
DESCRIPTION 

Random tree dbhs based on numbers of each dbh data segments in one hectare 
(data from Dr. Imhoff), and tree positions are generated with uniform dis- 
tribution in a stand, which is defined as an area of 10 hectares or 160 pix- 
els (each pixel with size 25 * 25 m). Radar returns are simulated for each 
model component in each pixel. The returns are further averaged to get 
returns of each pixel and of a stand. 

SIGMA is executable. 

intab file is a look-up table of radar cross sections of dielectric cylinders derived 
from exact solutions for Sundri and Gewa. This table covers one angle 
and tree dbh from 5.0 to 39.9 cm with 0.1 cm as a dbh step. 

indata file contains parameters needed for simulating a forest stand with optional 
ground conditions and optional solutions of model components. 

inpxl file is the input data, which has 14 columns of dbh segments (5 — 9.9, 10 
~ 14.9, 15 -- 19.9, 20 - 24.9, 25 - 29.9, 30 - 34.9, and 35 - 39.9) for 
Sundri and Gewa in one pixel corresponding to tree numbers of each dbh 
segment in each pixel. The tree positions or (x, y) coordinates are also 
randomly generated with uniform distribution in a stand. 

indbh file contains all the data (dbh' >= 5cm) with one data number per line for 
Sundri and Gewa in the whole stand. 

sts_o is a output file of statistical data (means and standard derivations) for each 
model component, and the sum of these components in a stand, respec- 
tively. 

t_o file summarizes radar returns of those components in each pixel. 

The rest *_o files are outputs corresponding to six radar backscattering com- 
ponents (1 to 6) mentioned above, respectively. 

NOTE 
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FILE pointers, *jptab, *fpin, *fppxl, *Jpdbh, *fps, *fpt, *fpl, *fp2, *fp3, *fp4, 
*fp5, and *fp6, point to intab, indata, inpxl, indbh, sts_o, t_o, dbs_o, 
vsu_o, vsl_o, tg_o, glc_o, and guc_o. 

The following are the subroutines called by the mainx program. 

head.c: defines all the external variables employed in the head.h file. 

indata.c: input parameters to choose model options and different surface conditions. Also, prints 
out these input parameters to check that the parameters are input correctly. 

openfl.c: open files to read and write. Totally, there are 12 files opened, in which four files are 
input and eight files are output. 

phs_x.c: random phases and/or x coordinates (distances in range direction) for Sundri and Gewa 
are generated by using randomO functions in the UNIX system, which are employed to count 
phase differences when summing each trunk / surface interaction terms. 

equi.c: calculate an equivalent value for two values x and y, weighted by their numbers nl and n2. 

ref.c: Fresnel reflection coefficient (both magnitude and phase) of H or V polarization incidence 
from eqns (9.1 - 49) and (9.1 -46) (Ruck, et al., 1970). 

gama.c*: incoherent component of forward scattering coefficient from a very rough surface using 
equations 9.7-118 - 132, in Ruck et al. v. 2, pp. 720-722. 

beta.c*: p parameters for a very rough surface scattering, which is defined in equations 9.1 - 124 
to 9.1 - 132 (p. 722, Ruck, et al., 1970). 

biomas.c: tree biomass in one pixel (25 * 25 m), which is calculated based on data from Dr. 
Imhoff. The biomass data is employed to compute attenuation coefficient of trees (Sundri 
and Gewa) for each pixel. Furthermore, the biomass data varies with pixels as a result of 

♦note: these programs are not used in this recent Bangladesh research project, but were written for the 

sake of completeness and for extension to different roughnesses of ground surfaces, in the future. 
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tree number change for each pixel. 

tau.c: the attenuation coefficient (t) of a canopy is calculated by using equations in Battan (1970, 
p. 67) eqns 6.5 and 6.6. The attenuation coefficient is equal to a scattering coefficient plus 
an absorption coefficient. 

ep.c: calculate dielectric constants (e’ and e") (note: e = e’ -j e") of mixture layers (the air as a host 
material, and materials, such as leaves, branches, trunks, etc. as inclusions) based on the for- 
mula in Ulaby, et al., (1986, E.62). The reflective model is chosen. 

vw.c: canopy volume scattering model (Attema and Ulaby, 1978). Further effort will be made in 
1989 to improve the simulation of the canopy component by extending this model component 
to include the disk model for broad-leaf evergreen forest (Eom and Fung, 1984; Lang and 
Saleh, 1985; Le Vine, et al., 1982, 1985). 

tge.c: radar forward scattering coefficient per unit area is derived from exact solutions of trunks 
(Sundri and Gewa), which are treated as dielectric cylinders with finite lengths and with 
smooth trunk surfaces. _ 

tgc.c: is the comer reflector model for tree trunks of Sundri and Gewa, which is the same as that 
used by Richards, et al. (1987). 

surl.c: backscattering coefficients per unit area from a slightly rough surface using a small pertur- 
bation model (Dobson and Ulaby, 1986), when the ground surface is in the water-off condi- 
tion. 

surll.c: backscattering coefficients per unit area from a very rough surface using eqs. 9.1-139 and 
141, in Ruck et al., vol. 2., p. 721 - 725. This program is again not used for the Bangladesh 
research project, but has been written for completeness in the whole modeling process. 

ee.c : interaction of canopy - surface scattering per unit area using Engheta and Elachi’s model 
(1982) by means of finding a° in terms of the water cloud model (Attema and Ulaby, 1978). 
This formula is similar to that used in Sun and Simonett (1988). 
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4. Accessory programs: 

The following are accessory programs, which are used to provide necessary input data in the 

correct format required by the above programs, and to maintain them handily. 

dbh.c: produces all the DBHs of total tree numbers (Sundri and Gewa) in a stand area, which are 
generated by the randomO function in the UNIX system with random uniform distribution. 
The input number is the total number of the dbhs*. 

pojree.c: random (x, y) positions of total Sundri and Gewa trees in a stand are generated by cut- 
ting the whole output data of the randomO function in the UNIX system into two parts. The 
first half of the data set is used to yield x coordinates, and the second half y coordinates. 
Twice the number of the positions are input. The (x, y) positions are designated in centime- 
ter units. The stand has an area of N hectares. 

po_num.c: uses the output of above (x, y) coordinate data as its input to produce the tree numbers 
of each dbh segment respectively for each simulated pixel in the whole stand. To run the 
program, tree number (TOTAL) for each dbh segment in a hectare and the position data 
should be in one input data file. The first line of input data is TOTAL and the rest of the 
data lines (N x TOTAL) are positions or (x, y) data, where N is a constant and can be 
changed for different stand sizes. In our case, the N is 10. 

makefile: maintains the above programs, which is based on the make command on the UNIX sys- 
tem. 

RUN: is an executable macro. Instead of typing the "SIGMA intab indata inpxl indbh sts_o t_o 
dbs_o vsu_o vsl_o tg_o glc_o guc_o" for each model being run, just type RUN, which does 
the same job as that inside the double quotation. 
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*note: all the dbhs range from 0.0 to 4.9 cm. Constants 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, and 35.0 are ad- 
ded to the dbh data set to produce dbhs in the following ranges (5.0 -- 9.9, 10.0 — 14.9, 15.0 — 19.9, 20.0 
— 24.9, 25.0 — 29.9, 30.0 — 34.9, and 35.0 — 39.9 cm) in the main.c program, respectively. 
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if no trees in a pixel, calculate the dbs component directly with 
the attenuation of the lower canopy only **/ 
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printf ("Sundrt biomass in a pixel: Ml - %lf\n“, Mis- 
print! (“Gena biomass in a pixel: M2 - tlf\n\n", M2); 
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if Si“-V‘ hh,th * t *' thetas ' phl ' epr - epi,; 

if P p Ph f_" 3 ^ thv(theta ' thetas, phi, epr, epi); 
btph2 - btvh (theta, thetas, phi, epr, epi); 
btph2 - btvv (theta, thetas, phi, epr, epi); 
- btph2 * ji; 
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printf ("dielectric constant :\n") ; 

print f ("epr - tlf, epi - -*lf\n\n", -epr, *epi); 
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for(l - 0; 1 < N_GEWA; ++P0, 1++) 

/* effective height of Gewa trunk */ 
ht - A3 + B3 * DBH_g [ 1 ] ; 
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